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Abstract 

We present a study of r~ — t- vr^vr+vr"!^,-, r~ — > J^^vr+vr^z^T-, t~ — t- K~ K'^tt^Ut, and 
T~ — )• K^K^K^Vr decays using a 666 fb^^ data sample collected with the Belle detector at 
the KEKB asymmetric-energy e^e~ collider at and near a center-of-mass energy of 10.58 GeV. 
The branching fractions are measured to be: B{t~ — )• Tr^Tr+Tr"/^,- )= (8.42 ± 0.00l°:|) X 10-2, 
B{t~ K~TT+TT~Ur)= (3.30 ± O.Oll[!;J^) X 10-3, B{t~ K~K+7r~Ur)= {1.55±OMtom) ^ 10~^ 
and B{t~ K~K+K~i^r)= (3.29 ± 0.17l:o;2o) x 10-^, where the first uncertainty is statistical 
and the second is systematic. These branching fractions do not include contributions from modes 
in which a vr+vr- pair originates from a decay. We also present the unfolded invariant mass 
distributions for these decays. 

PACS numbers: 13.35.Dx, 12.15.Hh, 14.60.Fg 
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INTRODUCTION 



The decays of the r lepton into three pseudoscalar particles can provide information on 
hadronic form factors, K* resonance spectroscopy, and the Wess-Zumino anomaly and 
also can be used for studies of CP violation in the leptonic sector By studying decays 
into final states that contain one or three kaons, it is possible to extract strange spectral 
functions, which can be used for a direct determination of the strange quark mass and 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \Vus\ js^Sl- These three hadron 
decays have been studied since the initial discovery of the r lepton, but only the decay to 
three charged pions has been closely investigated. Even for this mode, measurements of its 
branching fraction have only been provided by the CLEO jo] and BaBar j3] groups. For all 
other modes, there are still many unresolved problems. For example, the branching fraction 
for r~ — K^tx^tx^Ut decay recently measured by the BaBar collaboration [3] is significantly 
lower than the values from most previous measurements js]; there is large scatter in the 
measured central values, which is reflected in the large scale factor (2.1) applied by the 
Particle Data Group when forming their average js]. Intermediate two- and three-body 
resonance states in this decay have been previously studied but only with limited statistical 
precision, see, e.g., Refs jsHl. There is evidence for quite rich dynamics, with strong signals 
for the i^i(1270), -R'i(1400) and /('*(1410) resonances in the three-body final state; other 
resonances can be seen in various two-body sub-systems. Theoretically, two intermediate 
resonances, the p(770)'^ and i^*(892)'^ (or its excitations), are expected to contribute to r 
decay into K~ii+tx~u^: r" K~p{nQf{-^ tt+tt")//^ and r~ K*{%92f{-^ K-'K+)TX~Vr 



12|. These are shown in Figs[T]J^a) and (b), respectively. 

Here we present new measurements of the branching fractions for — )■ tt^tt+tt^z/^, 
—7- /i ^TT+TT^z/^, — !■ K~K'^tt~i't-, and — )■ K~ K^v^- decays. (Unless otherwise 
specified, the charge-conjugate decay is also implied throughout this paper.) Events in which 
a 7r"*"7r~ pair is consistent with a decay are excluded. Because of particle misidentifica- 
tion, measurements of all four of the above modes are correlated and, thus, are considered 
simultaneously. 

We use a data sample of 666 fb~^ collected with the Belle detector at the KEKB 
asymmetric-energy e^e~ collider 13| on the T(4S) resonance, 10.58 GeV, and 60 MeV 
below it (off-resonance). This sample contains 6.12 x 10® produced r-pairs. The Belle de- 
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FIG. 1: A schematic view of e~^e~ — t- r^r" in the center-of-mass system, where one r (r"*") 
decays to the signal mode (r^ — t- K'^ tt^ tt^ v-t) and the other r (t~) decays to a leptonic mode 
(r~ —7- e~Vri'e oi^ — ?• iJrvr^'y)- Each figure shows one of two possible intermediate resonant 
states in the signal mode. 



tector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located inside a su- 
perconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux return located 
outside the coil (KLM) is instrumented to detect K\ mesons and to identify muons. The 
detector is described in detail elsewhere [141 . 



SELECTION OF EVENTS 



Three-prong decays of the r are selected from candidate t^t~ pair events as follows. We 
use events where the number of charged tracks is four and the sum of their charges is zero; 
the transverse momentum of each track in the laboratory frame is required to be larger than 
0.1 GeV/c, and the tracks should extrapolate back to the interaction point within ±5 cm 
along the beam direction and ±1 cm in the transverse plane. The sum of the reconstructed 



5 



momenta in the center-of-mass (CM) frame is required to be less than 10 GeV/c, and the 
sum of energies deposited in the calorimeter is required to be less than 10 GeV. The largest 
of the four tracks' transverse momenta is required to be greater than 0.5 GeV/c to reject 
two-photon events, which typically have low transverse momentum tracks. To reject beam- 
related background, we require the position of the reconstructed event vertex to be less than 
3 cm from the interaction point along the beam direction, and separated by less than 0.5 
cm in the transverse plane. The missing mass, M^j^^ = (j^jnit ~ Xltr^'tr ~ Yli-yPi)"^^ ^^"^ 
polar angle of the missing momentum in the CM frame are efficient variables for rejecting 
two-photon and Bhabha backgrounds. In the definition of the missing mass, j>tr and are 
the four-momenta of measured tracks and photons, respectively, and j>init is the initial CM 
frame four-momentum of the e+e~ beams. We require that the missing mass be larger than 1 
GeV/c^ and less than 7 GeV/c^, and that the polar angle with respect to the beam direction 
in the CM frame be larger than 30° and less than 150°. Detailed information on the t^t~ 
pair selection and background suppression using missing mass can be found in 15|. 

Particle identification is used to select r events that contain one lepton (electron or muon) 
and three hadrons (pions or kaons). The magnitude of thrust is evaluated and is required 
to be larger than 0.9 to suppress two-photon and e^e" — )■ qq backgrounds, where the thrust 
is defined by the maximum of (^^ \n ■ 'Pi\)/(Yln Here pi is the momentum of the i-th 

track and n is the unit vector in the direction of the thrust axis - the direction maximizing 
the thrust. We require that the angle between the total momentum of the hadrons and the 
lepton momentum in the CM system be larger than 90°, since the tag-side lepton and signal- 
side hadrons usually lie in opposite hemispheres: the so-called 1-3 prong configuration, due 
to the large transverse momenta of r leptons. The invariant mass of charged tracks and 
gamma clusters on each side is required to be less than the r mass. Finally, we require that 
there be no K^, vr^, or energetic 7 on the signal side, where the selection criteria for these 
particles are described below. Figure [2] shows distributions of the variables used to select the 
1-3 prong sample, where in each case, the selection value is indicated by a vertical line in 
the figure. The surviving events are candidates for — hih^h^Vr-, where /ii,2,3 is either a 
pion or a kaon candidate. At this stage, the efficiencies for reconstructing — t- tt^tt+tt^z/^, 
r~ — )■ K^'k'^'k^Vt-, — K~ K^tt~Vt-, and r~ — K~ K~ decays as — t- h~h^h^Vr 
are 26.7%, 27.5%, 27.2%, and 24.8%, respectively, while the reconstruction efficiencies of 
the dominant background modes, — )■ 7r^7r+7r^7r'^z/^, r~ — )■ K^t\:~Vt-i and e^e~ — gg, are 
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6.0%, 1.8%, and 0.004%, respectively. The fractions of r~ — 5> 7r~7r+7r~7r''z/T-, r~ — )■ KgTr'Ur, 
e~^e~ — )■ qq and two-photon backgrounds that are reconstructed as — )■ h~h^h~Vr are 
6.47%, 0.34%, 0.35% and 0.05%, respectively; the two-photon background contribution is 
negligible. 




FIG. 2: Distributions of variables used in the event selection, (a) The angle between the total 
momentum of the hadronic system on the signal-side and the lepton momentum on the tag-side in 
the CM system, (b) The magnitude of the thrust, (c) The invariant mass of the hadronic system. 
The triangular points show data, the dashed histograms show the r~ — t- K~ti'^ti^Ut- signal MC, the 
solid histograms show the two-photon background, and the dotted histograms show the e'^e~ — >• qq 
background. The dotted vertical lines show the boundaries used to select the signal candidates. 
All samples are normalized to the same number of events. 

An important issue for this analysis is the separation of kaons and pions. In Belle, dE/dx 
information from the CDC, hit information from the ACC, and time-of-fiight measurements 
from the TOF system are used to construct a likelihood for the kaon (pion) hypothesis, 
L(A') (L(7r)). Figure [3] shows the kaon likelihood ratio, L{K) / {L{K) + L{n)), as a function 
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of momenta. Clean kaon-pion separation is evident over the entire momentum range relevant 
for this analysis, where the track momentum ranges from 0.1 GeV/c to ~5 GeV/c, and the 
average momentum is ~1.3 GeV/c. As discussed below, we choose a relatively stringent 
particle identification (PID) requirement for kaons {L{K)/{L{K) + L{tt)) > 0.9) and a less 
restrictive one for pions {L{K)/{L{K) + L^ir)) < 0.9). The kaon identification selection has 
an efficiency of ~73% for kaons and rejects ~95% of pions. The kaon and pion identification 
efficiencies are calibrated using kaon and pion tracks in data from kinematically identified 
D*^ decays, D*~^ — D^{^ A'"?!"^)??"*". We evaluate the efficiencies and fake rates for 
this calibration sample and compare them to Monte Carlo (MC) expectations. From this 
comparison, we obtain a correction table as a function of track momenta (piab) and polar 
angles {Oiat), and apply it to the MC. The accuracy of the correction factor, which is the 
source of the systematic uncertainty of the evaluation of the branching fraction and the mass 
spectra, is limited by the statistical uncertainties of kaon and pion samples from D*~^ decays 
in the certain piab and 6iab bins, and the uncertainty of the D*~^ signal extraction. 



FIG. 3: Likelihood ratio for kaon identification, L{K) / {L{K) + L(7r)), versus particle momentum. 
The filled and open circles correspond to kinematically identified kaons and pions, respectively. 

The kaon identification criterion is determined by maximizing a figure-of- merit (FOM), 
where the requirement on the kaon identification likelihood ratio is varied. The figure-of- 
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merit we used is defined as: 



S 



FOM = ^ , (1) 

where 5" is the number of signal events and is the number of cross-feed background 
events. For the r~ — )■ K'tt^ti'Vt signal (S*), the cross-feed background component (A^) 
includes contributions from r~ — t- -k^tx^ti^Vj. and r~ — t- K~ K^ti~Vt- decays. The result of 
FOM optimization for the — K^n^TT^Ur decay is shown in Fig. HJ where one can see 
that FOM is indeed maximal for the particle identification criteria used in this analysis. The 
same FOM analysis was performed for the r~ — tt^tt^tt^Uj. and r~ K~ K'^ti~Vj. decay 
modes and resulted in similar particle identification criteria. For electron identification, the 
likelihood variable is calculated based on the track extrapolation to the ECL, the ratio of 
the energy deposited in the ECL to the momentum measured in the CDC, the measured 
dE/dx in the CDC, the shower shape in the ECL, and light yield in the ACC 16|. The 
extrapolation of tracks from the CDC and SVD to the KLM is used to construct the likeli- 
hood variable for muon identification jl^]. The efficiencies and systematic uncertainties for 
lepton identification are evaluated by using a 77 — )■ e~^e~ / fi'^ fi' control sample. 
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FIG. 4: Figure-of-merit as a function of the likelihood ratio L{K) / {L{K) + L{tt)) requirement, used 
to separate the kaon and pion hypotheses. The uncertainties are evaluated from the comparison of 
the efficiencies and fake rates of kaon identification with those of the D*^ — t- D^{-^ K~tt^)'k^ con- 
trol sample. 



To reduce the feed-down background and to estimate the remaining backgrounds that 
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contain and tt*^ mesons, we reconstruct and vr'' signals explicitly. candidates are 
formed from oppositely charged pairs of pion tracks with invariant mass MijiTi) within ±13.5 
MeV/c^ (iSa) of the mass. To improve the purity, the point of closest approach 
to the interaction point along the extrapolation of each track is required to be larger than 
0.3 cm in the plane transverse to the beam direction. The azimuthal angle between the 
momentum vector and the decay vertex vector of the reconstructed is required to be 
less than 0.1 rad. The distance between the two daughter pion tracks at their point of 
closest approach is required to be less than 1.8 cm, and the flight length of the in the 
plane perpendicular to the beam direction is required to be larger than 0.08 cm. Finally, 
we require the invariant mass of each candidate to be inconsistent with that of a A, A, 
and a photon conversion, where the daughter tracks are assumed to be pions and protons 
(antiprotons) or electrons and positrons as appropriate. Candidate vr^'s are reconstructed 
from pairs of photon clusters. The energy of each photon is required to exceed 50 MeV for 
candidates in the barrel part of the calorimeter, and 100 MeV for candidates in the endcap 
part of the calorimeter. We also require the invariant mass of the two photons, M(77), to 
be within ± 16 MeV/c^ of the 11^ mass. Events containing one or more reconstructed or 
7r° are rejected. In addition, in order to reject events containing tt^'s, we require that there 
be no photon with energy larger than 0.3 GeV. 



EFFICIENCY ESTIMATION 



To estimate the signal efficiency, 97 million signal MC events for the r" — t- tt^tt+tt^z/t- 
decay mode and 5 million signal MC events for each of the — K^ir^-jT^Ur, — )■ 
K~ K^'K~jjj-, and — )■ K~K~^K~Ut- modes were generated using the TAUOLA 18| based 



KKMC 



GEANT3 



19 1 generator. The detector response for all MC data sets was simulated with the 
20| package. The TAUOLA generator provides models for r~ — )■ 7r~7r+7r~z/^, 



— 7- Ji ~7r"^7r~z/T-, and — )■ K~ K~^7r~Ur decays, where various intermediate resonances 
are taken into account. To check whether there is efficiency bias related to the specific decay 
model used, events with phase-space decay distributions were generated using the KKMC 
program. The efficiencies as functions of the invariant mass for both r~ — t- K^n^Ti^Ur decay 
models are compared in Fig. The relative difference in efficiency between the TAUOLA 
and phase-space decay models is around 1%. Similarly, very small discrepancies between 
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the TAUOLA and phase-space decay model are observed in the r~ — > 7r~7r"'"7r~z/^ and r~ — )■ 
K~ K^'k~Vt- decays. We used the TAUOLA decay model to evaluate the efficiencies and 
their dependencies on the invariant mass of — tt^tt+tt^z/^, — /i "tt+tt^z/t-, and r~ — 
K~ K^Ti~Vj. decays. For the r~ — K~ K~Vj^ decay mode, the TAUOLA generator does 
not provide any decay model, therefore the efficiency is calculated assuming a phase-space 
decay. In practice, the efficiencies are evaluated using a two step procedure. First, we 
evaluate the response matrix for unfolding the mass spectra using the TAUOLA decay 
model. Next, we recalibrate the efficiencies using the unfolded spectra. Details of this 
procedure are described in the next section. 



Efficiency (TAUOLA decay) 
Efficiency (pliase space decay) 




"5.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 

M(K7t7i) [GeV/e=] 



FIG. 5: The efficiency for the r K t:~^tt decay as a function of M{KTnr). Squares and 
triangular points represent the TAUOLA decay model and phase-space decay, respectively. 



The average efficiencies and the fractions of the cross-feed backgrounds for all three- 
prong decays are summarized in Tabled The probability of reconstructing r~ — )■ 7r~7r"'"7r~z/^ 
as — K^Ti^Ti^Ur Is relatively low, but due to the large branching fraction of — )■ 
TT^Ti^TT^Ur dccay, there is a substantial contamination of misidentified r~ — > tt^tt+tt^z/t- 
events in the r~ — )■ K'tt^tt'Vj. sample. For a similar reason, the r~ — K~ K^K~Vj. decay 
mode has a large cross-feed background from r~ — K~K^'K~Vr- 
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BRANCHING FRACTION AND MASS SPECTRUM CALCULATION 

The numbers of candidate events in the r~ — )> 7r~7r"'"7r~z/^, r~ — K~7r~^7i~i'r, t~ — )■ 
K~ K^'k~Vt-, and — )■ K~ K~Vt- modes, after applying all selection criteria, are sum- 
marized in Table [TTl Possible backgrounds in these signal candidate samples include (a) 
cross-feed from the signal modes and (b) other processes such as — Tr^Tr+TT^Tr^z/,-, 
T" —7- K^Ti^Vri as well as continuum e^e^ — t- qq. The fractions of backgrounds coming 
from other processes are 5% to 12%, as summarized in Table [TTl (fourth column) for each 
signal mode. The fifth column shows the fraction of the main background component as 
determined from MC. The mode r" — -k^-k^-k^'k^Vt- dominates for the r~ — t- tt^tt+tt^z/t- and 

— 7- K^-k^-k^Vt- modes, while the background from continuum e^e~ — )■ qq is dominant for 
the r~ — )■ K~ K'^'K^Vr and r~ — )■ K~K'^K~Vt modes. 

To take into account the cross-feed between the decay channels, the true yield Nf^^° {i 
is the mode number and the values from 1 to 4 correspond to — t- n^n^ir^Ur, — t- 
K^Ti^'K^Ur, r~ —7- K~ K'^TT^Ur, and r~ — K^K'^K^u-r), is obtained using the following 
equation: 

j 

TABLE I: Summary of the efficiencies and the fractions of cross-feed backgrounds. The values in 
parentheses are the differences of the efficiencies or fake rates from those evaluated directly from 



the MC with the TAUOLA decay model, {sff - in percent (see text). 





Generated decay mode 


Reconstructed 


T — )• 


T — 7- 


T — )• 


T — )■ 


decay mode 


TT^ TT^ TT^ ly-j- 






K-K+K-Vr 


—7- TT^TT^TT^Vt- 


0.22 (-0.1%) 


0.079 (-0.1%) 


0.022 (-1.1%) 


6.4 X 10-3 (+2.6%) 


T~ — >• K~'K'^T:~Vr 


0.012 (-0.2%) 


0.18 (-0.5%) 


0.047 (+1.4%) 


0.019 (-3.1%) 


T~ —7- K^K'^TT^Ur 


3.9 X lO-'' (-0.3%) 


4.7 X 10-3 (-4.4%) 


0.12 (-1.0%) 


0.051 (-0.4%) 




5.0 X 10"*^ (+4.5%) 


1.3 X 10-^ (-9.1%) 


2.3 X 10-3 (-14.0%) 


0.081 (+1.2%) 
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TABLE II: The number of reconstructed events (second column), the number of true events (third 
column), the number of background events other than three-prong cross-feeds (fourth column), 
and the main source of other backgrounds with its fraction of the total (fifth column). 



Decay mode 


^rcc 


^truc 




Main component in 






8.86 X 10^ 


3.52 X lO'^ 


9.35 X 10^ (10.6%) 


T~ — 


> vr vr^vr tt^Ut 


(64.2%) 


T~ — )• K~'K^TT~Vr 


7.94 X 10^ 


1.38 X 10^ 


9.60 X lO'' (12.1%) 


T~ — 




(34.4%) 




1.08 X 10^ 


6.47 X 10^ 


7.16 X 10^ (6.66%) 


e~^e~ 


qq (30.3%) 




T- K-K+K-i^r 


3.16 X 10^ 


1.37 X 10"^ 


1.71 X 10^ (5.41%) 




^ qq (53.0%) 





where Nj°^ is the number of events for the j-th reconstructed decay mode and N° is the 
number of background events in the j-th mode. Here Sij is the efficiency for detecting the 
j-th mode as the z-th one. The values of A^'''^'^ and N°^^'^'^ are hsted in the second and third 
columns of Table [Tll respectively. 

We have determined the efficiency (migration) matrix iteratively. First, it is determined 
using the mass spectra as modeled in the current TAUOLA MC program. (We refer to this 
efficiency matrix as ^Ij^)- As is discussed in the next paragraph, we then unfolded the mass 
spectra of each mode to obtain the true mass distribution. The resultant unfolded mass 

(2) 

spectra are then used to determine the efficiency matrix S!>j using a weighting procedure. 
The matrix elements S^j^ are given in Table [H where the values given in parentheses are the 
percentage differences between sj;J^ and Slj^; [S^j'^ —£'lf )/£^\ Since the efficiency is smooth 
as a function of hadronic mass, the difference between E^^^ and E^J^ is very small (~1%) in 
most cases, except for the off-diagonal entries related to the r~ — K~ K~Vt- mode. The 
source of these large discrepancies is the mass dependence of the fake rates, together with 
the considerable differences between the unfolded spectrum and the generated spectrum. 

The unfolding procedure extracts the true mass spectra for all four modes, denoted by 
the vector x, by solving the matrix equation Ax = b. In this equation, b is the vector of 
the reconstructed mass spectrum with backgrounds other than the three-prong cross-feed 
subtracted, and A is the response matrix. 

The matrix A takes into account the cross-feed between different modes caused by particle 
misidentification, as well as the effects of finite resolution and the limited acceptance of the 
detector. The response matrix is obtained by merging 16 sub-response matrices as shown in 
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Fig. |6l The matrix is determined using the TAUOLA decay model MC simulation. Clear 
correlations between the measured and generated values are seen along the diagonal parts 
of i. 
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FIG. 6: Response matrix A obtained by merging 16 sub-response matrices that represent the 
efficiency or fake rates among the four three-prong modes. The bin numbers correspond to the 
fohowing ranges of the invariant masses: (1) from to 138 , 0.41 GeV/c^ to 1.8 GeV/c^ of M(7r7r7r); 
(2) from 139 to 241 , 0.77 GeV/c^ to 1.8 GeV/c^ of MiKmr); (3) from 242 to 309 , 1.12 GeV/cP 
to 1.8 GeV/c^ of M{KKtt); and (4) from 310 to 341 , 1.48 GeV/c^ to 1.8 GeV/c^ of M{KKK). 
Ah bin sizes are 10 MeV/c^. 



For the unfolding algorithm, we follow the method of the ALEPH Colla 
based on the Singular Value Decomposition technique for matrix inversion 



Doration, which is 
2l| and advanced 
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regularization 22|. Two independent r~ — 5> h~h^h~Vr signal MC samples are used to check 
the validity of the unfolding procedure. One sample is used for the determination of the 
response matrix. We check the reproducibility of the generated M{hhh) mass spectra with 
the other sample, and find, for the whole mass range, that the differences are less than 1%. 
These discrepancies are included in the systematic uncertainties of unfolding (UNFl below). 
The resulting unfolded normalized mass spectra of the three-prong decays are shown in 
Fig. [71 where the error bars indicate the statistical uncertainty only, and the gray bands 
correspond to the systematic uncertainty. Note that the uncertainties common for all bins 
are taken into account as the systematic uncertainties on the branching fraction (Table HVT) . 
while the bin-by-bin errors remained in the normalized mass spectra l/N{dN/dM) are 
accounted here. The following sources of systematic uncertainties in the mass spectra are 
considered: the unfolding procedure (UNFl and UNF2), the kaon identification efficiency 
(KID), the background estimation other than the three-prong cross-feed (BGE), the effect 
of 7 veto (GAM), and the track momentum scale (MOM). Table UTTl summarizes various 
contributions to the systematic uncertainties on the normalized mass spectra l/N{dN/dM) 
for each decay mode. 

The systematic uncertainty due to the unfolding procedure is determined from MC by 
comparing the generated and unfolded mass spectra (UNFl). Another estimation of the 
uncertainty of the unfolding procedure is obtained by changing the value of unfolding pa- 
rameter that determines the effective rank of the response matrix (UNF2). This is the most 
important systematic uncertainty for KKK mode. The kaon identification efficiency and 
the fake rate are calibrated by using D*'^ sample as described in the previous section. The 
uncertainties due to the particle identification (KID) are evaluated by varying efficiency 
and cross-feed fraction by one standard deviation in the response matrix. The systematic 
uncertainty from the background estimation (BGE) is evaluated by varying the branching 
fraction values used in the MC by one standard deviation. The uncertainty due to the 7 veto 
(GAM) is evaluated by using different selection criteria for photon energy and comparing 
the resultant unfolded spectra. For the uncertainty from the track momentum scale (MOM), 
we reconstruct (f) — j- K^K~ and -D^ — )■ /T^tt^tt^ decays. By comparing the reconstructed 
masses of and D"^ mesons with their world average values, we conservatively assign 0.01% 
uncertainty for the momentum scale. By applying this variation to the reconstructed mass 
spectra of r~ — h~h~^h~h'T- decay and unfolding it, we obtain the uncertainties due to the 
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momentum scale. 

The systematic uncertainty varies as a function of the mass of the hadronic system in 
each mode. The typical uncertainties averaged over all mass bins are evaluated by summing 
the systematic uncertainties of each bin and ignoring the correlations between bins, and 
the results are 0.7%, 2.2%, 2.2%, and 9.5%, for the r~ tt^tt+tt^z/^, A'^tt+tt^z/^, 
r~ — )■ K~ K^TT~Vri and K~ K" Vr decay modes, respectively. In practice, we have 

obtained the covariance matrix for each source of the systematic uncertainties and added 
them to obtain the total uncertainty. The gray bands in Fig. [7] show the square roots 
of the diagonal components of the covariance matrix. The off-diagonal components of the 
covariance matrix are shown in Fig. [S]in term of the correlation coefficients p^^'^'^{i,j) defined 
by p^P^'^(i,j) = cov(z, j)/ ^Jcaw{i,i) ■ cov(j,j), where i and j are bin numbers. 



Figure [7] shows rather large discrepancies between the mass spectra of the genera- 
tor and our unfolded spectra except for r~ — )■ tt^tt^tt^z/t- decay. For example, for the 
r~ — K^n^n^Ur mass spectrum (Fig. [7b), it is evident that there are contributions from 
two resonant peaks, similar to the MC expectations. However, the faster rise of the real 
data as compared with the MC in the 1.0 — 1.1 GeV mass region and the fall-off at high 
masses in the r~ — A'~7r+7r~z/^ mass spectrum may also be suggestive of a substantial 
non-resonant component that is not included in the current MC model. The use of the 
unfolded invariant mass distributions in estimating the efficiencies allows one to compensate 
for model deficiencies when evaluating the branching fractions. 

In order to determine the branching fraction, we normalize to the number of pure leptonic 



decays, where one r decays to r — )■ eVv and the other decays to r — )■ fiVu j23| (hereafter such 
events are referred to as { e, /i } events). The branching fraction for the i-th decay mode can 
then be written as: 

where N^i^^^,^ = 8.1 x 10^ and £e^j = 0.221 ± 0.003 are the number of {e,/i} events and 
the corresponding detection efficiency, respectively. The uncertainty in eg^i is determined 
from the lepton identification uncertainties. Here Br^evv = 0.178 and Br^^^vu = 0.174 
are the branching fractions of r — eUu and r -> pViy decays, respectively, and Br^wv = 
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TABLE III: Summary of the relative errors of the unfolded mass spectra, 1 /N{dN/dM) (in %) from 
different sources of uncertainties: the unfolding procedure (UNFl, UNF2), the kaon identification 
(KID), the background estimation (BGE), the 7 veto (GAM), and the momentum scale (MOM). 
The "average" uncertainties (the 2nd, 4-th, 6-th, and 8-th columns) are evaluated by taking average 
of errors in all bins. The "peak" uncertainties (the 3rd, 5-th, 7-th, and 9-th columns) represent 
the errors at the peak position of the spectra. See the text for a more detailed description. 



Sources 


T — 




T — 


> 


r — 




r — 




of 


TTTTTTU 


Kttttu 


KKttv 


KKKu 


uncertainties 


average 


peak 


average 


peak 


average 


peak 


average 


peak 


/'TTATT?1 \ 

( U IN r' i j 


0.5 


0.5 


0.1 


0.1 


0.4 


0.4 


0.6 


0.6 


(UNF2) 


0.1 


0.1 


1.6 


1.5 


0.9 


0.9 


7.3 


4.0 


(KID) 


0.4 


0.4 


1.0 


0.8 


1.5 


1.1 


1.9 


0.9 


(BGE) 


0.2 


0.2 


0.6 


0.5 


0.4 


0.3 


1.9 


1.1 


(GAM) 


0.6 


0.4 


0.9 


0.4 


1.2 


0.9 


4.6 


3.2 


(MOM) 


0.1 


0.0 


0.4 


0.2 


0.3 


0.3 


2.7 


3.1 


Total systematics 


0.7 


0.7 


2.2 


1.9 


2.2 


1.7 


9.5 


6.2 


Statistics 


0.2 


0.2 


1.1 


0.8 


1.0 


0.8 


7.1 


4.6 


Total 


0.9 


0.8 


2.4 


2.0 


2.4 


1.9 


11.8 


7.7 



ISr^evu + ISr^^vu- The tiuc yields N^^'^'^ obtained from Eq. [2] are given in the third column 
of Table M 

This normalization method requires a precise measurement of the number of { e, /i } events 
and a careful determination of the corresponding efficiency. Figure [9] shows a comparison 
of the invariant mass of the electron and muon system and the cosine of the angle between 
electron and muon, for real { e, } events and the sum of MC expectations. The MC 
reproduces the data reasonably well. The background fraction is estimated to be 5.9% using 
MC. After subtracting background, we obtain the number of { e, /i } events. 

Figure [10] shows the observed distributions of M(7r7r7r), M{K7T7t), M{KK7t), and 
M[KKK). The backgrounds are represented by the hatched and the gray histograms for the 
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FIG. 7: The unfolded mass spectra of the three-prong decays, l/N{dN/dM): (a) M^mnr) distribu- 
tion of —7- Tr~^ TT^ I'r , (b) M{K'K'k) distribution of t" — )■ K~ ■k'^ tt~ , (c) M{KKir) distribution 
of r~ —7- K~ K'^TT^u-T, and (d) M{KKK) distribution of r~ — )■ K'^ K~Ut-. The black points 
correspond to the unfolded mass spectra with statistical uncertainties only, and the gray bands cor- 
respond to the systematic uncertainty. The solid histograms are the model predictions implemented 
in the current TAUOLA MC simulation. 

sum of the three-prong cross-feed backgrounds and other background components, respec- 
tively. Note that the three-prong cross-feed backgrounds are determined by the unfolding 
procedure previously described. 

Table HVl summarizes various contributions to the systematic uncertainties on the branch- 
ing fractions. The uncertainty in the track-finding efficiency is estimated from a comparison 
of real and MC data for the D* 7rD°, D° ■KTxK^g{K^g tt+tt^) decay sample. The 
uncertainties due to the kaon and pion identification are evaluated by applying the un- 
certainties in kaon identification efficiencies to the efficiency matrix in Table [B First, we 
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FIG. 8: The correlation coefficients of the unfolded mass spectra for three-prong decays. See 
Fig. [6] for the convention of the bin numbers. The pattern is the reflection of the detailed unfolding 
procedure. 



measure the kaon and pion identification efficiencies and their uncertainties using the con- 
trol sample of D*^ — )■ -D°vr+ and — )■ K~tt~^ events. Since these efficiencies are applied 
to the efficiency matrix, we propagate the uncertainties in the kaon and pion identification 
efficiencies to the covariance matrix of branching fraction measurement. The kaon identifi- 
cation systematic uncertainties for r~ — Ti^ii^Ti^Ur, — )■ K^ir^ir^Ur, r~ — K'^ K~^tt~i't-, 
and — )• K~ K'v^ decays are 1.3%, 3.9%, 1.9%, and 5.4%, respectively. For lepton 
identification efficiency uncertainties, we use 77 — )■ e~ / fi~^ fi~ events. 

In order to take into account the effects of the mass spectra on the branching fraction 
measurements, we have determined the mass spectra of the four dominant decay modes 
iteratively. The uncertainties in branching fractions due to the mass spectra are evaluated 
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FIG. 9: Distributions of the {e,fi} events: (a) Invariant mass of the electron-muon system, (b) 
Cosine of the angle between the electron and muon. The closed circles and the solid histogram rep- 
resent the data and the sum of the MC expectations, respectively. The open histogram represents 
the { e, /i } signal events, while the r-pair background and the two-photon background are shown 
by the hatched and the gray histograms, respectively. 



by varying the unfolded spectra by one standard deviation of their statistical and systematic 
errors, assuming 100% bin-by-bin correlations for systematic errors. The trigger efficiency 
for the three-prong modes is ~86% with an uncertainty of ~0.6%. This efficiency is obtained 
using a trigger simulation program applied to the signal MC samples. The uncertainty due to 
the 7 veto is evaluated by using different selection criteria for photon energy. The uncertainty 
due to the Kg veto is evaluated by removing the veto requirement. The uncertainties in 
background estimation include the effect of luminosity, the cross section for e~^e~ — r+r~, 
and the subtraction of background from non-three-prong r decays and continuum e~^e~ — )■ 
qq. The uncertainty on the r decay background subtraction is evaluated by propagating 
the errors from the r branching fractions other than those with three-prongs. The fraction 
of the e^e~ — t- qq background is negligible for — Ti^n^n^Ur and — K^ir^ir'Ur 
modes, and is 2.0 ± 0.2% and 2.9 ± 0.3% for r" ^ R-R+Ti-Vr and K~K^K-Ur 
modes, respectively. The uncertainties in the e^e~ qq background are estimated from 
comparison of the number of events in data and e^e~ — ?■ qq MC, in the mass region above 
the r mass. The effects of the uncertainties in the luminosity and the cross section for 



24j are negligible because we use { e, yU. } events for normalization. The 



uncertainties in the leptonic decay branching fractions for the r are also taken into account. 
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FIG. 10: The reconstructed invariant mass distributions: (a) M(7r7r7r) distribution for t~ — )■ 
7r~7r+7r~z/^, (b) M^Kmr) distribution for — )• iir~7r~^7r~i/^, (c) M{KKtt) distribution for r~ — )• 
K'^TT^i'r, and (d) M{KKK) distribution for t~ — )• K~K'^K~Ut-. For all histograms, the black 
points with error bars represent the data, the gray histograms are the sum of all backgrounds other 
than the three-prong cross-feed, and the hatched histograms are the sum of three-prong cross-feed 
backgrounds. Note that the cross-feed backgrounds are estimated from the unfolding of the data. 
The uncertainties shown for the data are statistical only. 



After taking into account the backgrounds, the efficiencies and various sources of sys- 
tematic uncertainties discussed above, we obtain the following branching fractions for three- 
prong decays: 
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TABLE IV: Summary of the systematic uncertainties on the branching fractions (%). 





r — > 


r — > 


r — 7- 


r — > 


Source 




Kttttv 


KKttu 


KKKu 


Tracking efficiency 


+2.2/-2.0 


±2.1/-2.0 


±2.1/-2.0 


±2.1/-1.9 


Particle ID 


±1.9 


±4.0/-4.1 


±2.3 


±5.4/-5.6 


Mass spectrum 


±0.1 


±0.1 


±0.1 


±0.5/-0.8 


Trigger efficiency 


±0.5 


±0.5 


±0.6 


±0.6 


7 veto 


±0.9 


±1.5 


±1.7 


±0.5 


K'l veto 


±0.2 


±0.2 


±0.1 


±0.3 


Background estimation 


±0.3 


±1.3 


±0.3 


±0.4 


Leptonic branching fraction 


±0.2 


±0.2 


±0.2 


±0.2 


Total 


±3.1/-3.0 


±5.0 


±3.6 


±5.9/-6.0 







= (8.42 ± 0.00(stat. 


)^°;|(sys.))xlO-^^ 


B{t~ — K^Ti^Ti' 




= (3.30 ± 0.01 (stat. 


)^°:i?(sys.)) X 10-3 , 


B{t- K~K+n- 




= (1.55±0.01(stat. 


)t'omi^YS.)) X 10-3 , and 


B{t^ K-K+K- 




= (3.29±0.17(stat. 


)l°i9(sys.)) X 10-5 . 



The covariance matrices for the branching fraction measurements are given in Table IVl in 
terms of the correlation coefficients p^^{i,j) = cov{i,j)/^ycov(i,i) ■cov(j,j), where i and 
j correspond to one of the four three-prong decay modes. There is a strong correlation 
between r" — t- TT^ir^TT^Ur and r" K^ir^n^Ur decays, which mostly comes from the effect 
of kaon identification. 

Using the results of the branching fractions and the correlation matrix, the ratios of the 
branching fractions involving kaons to the branching fraction of — ?■ ir^TX^'n'^i'T- decay are 

^ The actual value of the statistical uncertainty of B{t^ — > tt^tt+tt^z^t-) is 0.003 x 10^^. 
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TABLE V: The correlation coefficients in the branching fraction measurements. 





T — 7- 


r — 7- 

K-K+TT-Vr 


r 

K-K+K^Ur 


T~ 


TT 




+0.175 


+0.049 


-0.053 


T~ 


K 






+0.080 


+0.035 


r~ 


K 








-0.008 



estimated: 

B{t~ K~7i^n~Ur)/B{T~ Ti-Ti+Ti-u^) = (3.92 ± 0.02+°;^^) x 10"^ , 

B{t^ K^K^TX^Ur)/B{T^ n-n^n~iy^) = (1.84 ± 0.01 ± 0.05) x 10"^ , and 

B{r- K~K+K-Vr)/B{T^ tt-tt+tt-z/,) = (3.90 ± 0.20l°:|) x 10"^ , 

where the first and the second uncertainties are statistical and systematic, respectively. 
DISCUSSION 

The results of this analysis for the branching fractions of various three-prong modes are 
listed in Table IVll together with recent results from BaBar [3]. Note that for the r~ — )■ 
TX^-K^Ti^Vj. and r~ K^tt^tt^Ut modes, the branching fractions listed do not include any 
contribution. 



TABLE VI: Comparison of the branching fraction results 



Decay mode 


BaBar 


Belle 


T~ — TT~ TT^ -ir^ Ur , % 
T~ — K~TT^TT~h>r, % 
T- K'K+TT-Ur, % 


8.83 + 0.01 + 0.13 
0.273 + 0.002 + 0.009 
0.1346 + 0.0010 + 0.0036 
1.58 + 0.13 + 0.12 


8.42 + o.oollj:^^ 

0.330 + 0.00llj];[JJ^ 
0.155 + 0.00ll[j:[jg^ 



In Fig. [m our results are compared with the previous measurements. For all modes 
studied, the precision of the branching fractions for both BaBar and Belle is significantly 

23 



higher than before. The accuracy of our results is comparable to that of BaBar, but the 
central values show striking differences in all channels other than — ?■ tt^tt+tt^z/t-. For the 
r~ — tt^tt+tt^z/t- mode, our result is lAa lower than that of BaBar, while for other modes 
the branching fractions obtained by Belle are higher by S.Ocr, 3. Oct, and 5Aa than those 
of BaBar for the — K^ti^tt^u^, — t- K~ K^ti^Vt, and r~ — K~ K'^ Vj. modes, 
respectively. 

CONCLUSION 

Using a data sample of 6.12 x 10* r^r^ pairs collected with the Belle detector, we measure 
the branching fractions for the — tt^tt+tt^z/^, — t- A' ^tt+tt^z/^, — )■ K~ K^Ti~Vri and 

— K~ K^K~Vj. decay modes. We have been able to extract unfolded invariant mass 
spectra for all four modes (Fig. [Tj) by taking into account the cross-feed effects. A future 
detailed analysis of these spectra as well as those of the intermediate two-body states will 
allow studies of the decay dynamics to be performed. 
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FIG. 11: Summary of the branching fraction measurements of three-prong r decays. 
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